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Abstract 

A D5 brane winding around a stack of D3 branes can be used as a model of persistent 
current in a thin superconducting ring, with the number N of D3s corresponding to 
the number of transverse channels in the ring. We consider, in the large N limit, 
classical and quantum stability of a superconducting state with a uniform winding 
number density q. We find that there is a classical superconducting solution with any 
q, no matter how large, but for q larger than a certain q^ it is only metastable and 
can decay into a state with a smaller q by quantum phase slips. We argue that the 
decay produces strings via a version of the Hanany-Witten effect (in a non-transverse, 
non-supersymmetric arrangement of branes). This parallels the requirement of quasi- 
particle production in a clean (disorder-free) wire in field theory of superconductivity. 
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Introduction 





In this paper, we consider solutions of type JIB supergravity that describe one brane 
winding around another. The specific example we choose corresponds to the configu- 
ration in string theory that consists of a stack of D3 branes and a single D5 brane, 
in the following arrangement: 

0123456789 
N D3s: X X X X (1) 
D5: XX X X X X 

where crosses denote the coordinates wrapped by the branes. The only spatial direc- 
tion shared by all the branes is x = x^, which we assume to be a circle, possibly of a 
large radius. We will be interested in the case when, as one goes around the circle, 
the D5 winds around the D3s, the total of W times. 

Our interest in this configuration is that we believe it to be suitable for description 
of a current- carrying state in a thin superconducting ring, with W corresponding to 
the winding number of the order parameter. There are questions concerning thin 
superconductors, such as, for instance, conditions of quantum stability of the super- 
current, that are not easily answered by the conventional mean-field theory, and we 
would like to see if supergravity can help us to answer these questions. 

At small values of the 't Hooft coupling A = QsN (where Qs is the closed string 
coupling), one may start by considering the trivial embedding — with the D3s at x^ = 
. . . = x^ = 0, and the D5 at x"^ = x^ = x^ = x^ = — and developing a field theory 
around it. This corresponds to replacing strings connecting the D5 to the D3s with 
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their ground states and results in a non-supersymmetric gauge theory of species 
of left- and right-moving massless defect fermions. These correspond to N channels 
of conduction electrons in a wire|l| 

The role of interactions, on the other hand, is most readily understood in the 
opposite, strongly coupled limit, when both and A are large. In this case, it is 
possible to replace the D3 branes by their classical geometry [T] while considering the 
D5 as a probe. The near-horizon region of the D3 geometry has a field-theoretical 
interpretation even at strong coupling: it is dual to an SU{N) superconformal gauge 
theory [2], [3[ H]. The trivial embedding of the D5 crosses the horizon but, as we will 
see, is unstable: once displaced, the D5 is pushed out of the near-horizon region, to 
large absolute values of 

+ ix^ = Ae^<^ . (2) 

The correspondence of [2|, [3l H] is not applicable at large A, so supergravity is the 
only available microscopic description of this strongly coupled system. It will be our 
starting point in what follows. 

We interpret the instability of the trivial embedding as formation of a supercon- 
ducting state [5]. The latter is represented by a nontrivial embedding, which we find 
numerically for various values of the winding of the phase 0. A nonzero winding cor- 
responds to a nonzero supercurrent in the wire. All embedding solutions we consider 
here are uniform in the x direction and so describe a superconductor in the clean 
limit (i.e., without disorder). In particular, the gradient of (f) (the winding number 
density, in units of l/27r) is independent of x and given by 

g = ^, (3) 

where L is the length of the wire. In practice, disorder is often important, and the 
first step towards including it in our description would be to consider x-dependent 
background geometries. We discuss this further in the concluding section. 

We find that there is a nontrivial embedding solution for any value of g, no matter 
how large. For each of these solutions, the minimal distance from the D5 to the D3s — 
the quasiparticle gap in units of the string tension — is nonzero. This is in contrast to 
the case of fixed current (and zero winding) we have considered in [5]. In that case, 
the superconducting solution is gapless at any nonzero current and disappears at a 
critical current. The distinction between the fixed- winding and fixed-current cases 
has a parallel in the conventional mean-field theory, as we discuss in Sec. [2l 



t In a wire, different channels correspond to different transverse wavefunctions and different 
projections of spin. Tfius, at large N , the material we are describing is one-dimensional only in its 
superconducting properties (e.g., the value of the order parameter), while the electrons are free to 
move in all three directions. 



3 



The unconditional classical stability of the supercurrent may seem surprising: the 
prevalent situation in superconductivity is when the supercurrent is classically stable 
only up to a certain g = g^, at which it becomes energetically favorable to produce 
quasiparticles with momenta antiparallel to the flow. The stability bound is known 
as Landau's criterion; we recall its derivation in Sec. [2l Note, however, that our 
superconductor is perfectly uniform and therefore momentum- conserving, while the 
Landau process is not. So, it is possible that the classical instabihty will not show up 
until we break momentum conservation by going over to an x-dependent background. 

Next, we consider quantum stability of the supercurrent, in particular, the role 
of quantum phase slips (QPSs). A QPS is a tunneling event that removes a unit of 
winding [S]. Since that releases momentum from the supercurrent, a suitable sink of 
momentum is required for the process to take place. Unlike in the Landau process, 
however, that sink may be a part of the system itself, so that no momentum needs to 
be supplied from the outside (e.g., by disorder). The existing microscopic theory of 
this effect [7] is based on the Bogoliubov-de Gennes (BdG) equations, which, when 
interpreted hberally enough (as reviewed in Sec. H]), allow one to go one step beyond 
the conventional mean field, to include fluctuations of the order parameter. The result 
is that, in a uniform system, a QPS must produce a certain number of fermionic 
quasiparticles, and the associated energy cost may block phase slips altogether. 

The liberal interpretation of the BdG equations, referred to above, treats the 
order parameter and the quasiparticles as separate quantum fields. One may be 
concerned that this involves double counting, as both are ultimately made of the 
same electrons. One may counter this by arguing that, as in any two-fluid model, the 
two fields simply represent two different transport mechanisms (in our case, of the 
electric charge), and the separation of the electrons into superconducting and normal 
should not be construed too literally. Nevertheless, it is of interest to develop an 
approach in which the superconducting and normal components are represented by 
different types of excitations — e.g., an elementary field vs. a soliton or a string. Dual 
gravity provides precisely this type of approach. 

An important role in our calculation is played by the Wess-Zumino (WZ) term in 
the D5 action. We find that for any nonzero winding W it causes the D5 to have an 
amount, equal to NW, of the charge conjugate to the worldvolume gauge field. Since 
a QPS changes W, it changes the amount of charge and, by charge conservation, the 
difference has to be picked up by excitations. We interpret this transfer of charge as 
a version of the Hanany-Witten (HW) effect [Sj: when a D5 brane crosses a D3, a 
fundamental string is produced. This version of the effect is similar but not identical to 
the one well known in the literature, in which the D3 and D5 have no common spatial 
directions and preserve some supersymmetry [9l [TOl [11]. We find it remarkable that 
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the present version completely parallels the requirement of quasiparticle production 
obtained in field theory of superconductivity on the basis of the BdG equations. 

While a QPS reduces the energy of the supercurrent, to see if it is in fact energet- 
ically allowed, we must also take into account the energy of the produced quasiparti- 
cles. We find that decay of supercurrent into quasiparticles via a QPS is energetically 
forbidden until the winding number density reaches a certain q^. Thus, although a 
solution with nonzero supercurrent exists for any g, at g > it is only metastable. 
We discuss some consequences of this in the concluding section. 

2 A review of the BdG formalism 

The BdG formalism is not used for any calculations in this paper. We review it here 
as a reference point for later comparison with the results obtained using gravity and 
also to illustrate the distinction between the fixed-winding and fixed-current cases. 

Consider the following Lagrangian, describing species of (1+1) dimensional 
fermions, ipn) which interact with an order parameter field $: 



Here a = 0, 1 and n = 1, . . . ,N . The superscript c denotes charge conjugation, and 
/X is a chemical potential. Each ipn is a two-component Dirac spinor. We use the 
representation of the 7 matrices with 7° = ai, 7^ = — io"2, where as are the Pauli 
matrices. Thus, in the normal state ($ = 0), the upper components of ipn describe 
right-, and the lower left-, moving fermions. 

The Dirac fermions ifjn can be related to the electron operators a^^L of a super- 
conductor by taking even and setting, for each v = 1, . . . , N/2, 



where an arrow denotes the projection of spin. The minus sign in the second in- 
stance in ([5]) causes $ to couple to a spin singlet, the case for most conventional 
superconductors. The index v labels the transverse channel. It refiects the fact that, 
while the fields ipn are functions of time and x only, the electrons in a superconductor 
move in all three spatial directions and can be classified according to their transverse 
wavefunctions. 

The BdG equations are obtained from (jlj) by variation with respect to ipn- Our 
interpretation of these equations will be broader than the conventional one: in most 
applications $ is considered as a classical background, but we will allow quantum 
changes in $. This is necessary if we wish to use (jl]) to describe quantum phase slips. 
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In mean-field theory, one can describe a steady supercurrent as a state with a 
uniformly wound order parameter: 

$(x) = <l>oe*«% (6) 

where $0 7^ is a constant. A quantum phase shp is described by a configuration 
(instanton) that interpolates in the Euclidean time between the state and a similar 
state with a smaller q, corresponding to one fewer unit of winding. Such an instanton 
has one zero mode for each fermionic species [121 [13]. An adaptation of the argument 
of [H] then shows that fermions must be produced in the phase slip process [7|. 

The quasiparticle spectrum near the ground state can be conveniently obtained 
by redefining the Fermi fields as follows: 

^e^'?"V^„(x,t). (7) 

This unwinds $ into a constant, ^(x) — > $0; but also modifies the Lagrangian: 

Lf ^ Lf - -J2 ^i^3^n ■ (8) 
^ n 

The additional term in (|8|) gives different chemical potentials to the right- and left- 
moving fermions. The corresponding Hamiltonian can be diagonalized by a Bogoli- 
ubov transformation. The resulting excitation energies are 

e^{k) = [{k-f,r + \^o\V±l- (9) 

The upper branch, e+(fc), describes excitations with momentum k, and the lower 
branch, e_(/c), those with momentum —k. Both reach minima a.t k = fi, which can 
therefore be identified with the Fermi momentum: ix = kp- The lower branch touches 
zero when q becomes equal to 

gL = 2|$o|, (10) 

that is twice the quasiparticle gap. Once q exceeds gL, it becomes energetically favor- 
able to produce excitations with momenta near —kp, and the ground state becomes 
unstable. The stability condition g < is known as Landau's criterion. As we have 
noted in the introduction (and discuss further in the conclusion), we do not recover 
this condition in the gravity-based calculation, presumably because of the perfectly 
momentum-conserving nature of our system. 

We can now define two limiting types of current- carrying states. In one type, 
quasiparticles are absent but the order parameter is wound as in ([6|); the current this 
state carries is a supercurrent. We refer to this C3iSG clS fixed winding. In the other 
type of state, (^{x) = $0, but the upper branch of the excitation spectrum Q is 
filled with quasiparticles, up to a certain finite density of them. The current these 
carry is a normal current. Because the quasiparticles form a Fermi surface ji| this state 

■''Not to be confused with the original electron Fermi surface at fc = /i. 
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is gapless. We identify it as a mean-field counterpart of the gapless superconductor 
described via a gravity dual in [S] . We refer to this CclSG clS fixed current. 

In general, one must allow for both supercurrent and normal current components 
to be present. As long as phase slips are neglected, the momenta of these components 
are separate conserved quantities. Phase slips, however, convert winding of the order 
parameter into momentum of quasiparticles and vice versa. The equilibrium winding 
and density of quasiparticles are then determined by comparing the free energies of 
different states with the same total momentum. In this paper, we will be interested 
in a condition under which the supercurrent-only state becomes metastable, and its 
decay begins to populate the normal component. 

3 Computation of the WZ term 

Equations of motion for the D3 geometry can be derived from the action 

s,, = J-Ji-.V^(n-^Gl,). (11) 

Here TZ is the Ricci scalar, and (7(5) is the Ramond-Ramond 5-form field strength. 
We use the shorthand notation 

ri2 _n n ABODE o^ 

*-^(5) ~ '^^ABODE'^ ■ U^j 

(^(5) is self-dual, which implies G^g^ = but, as is common practice, we impose 
self-duality in the equations of motion, rather than directly in the action. 

A stack of N extremal D3 branes is a solution to the equations of motion following 
from ( ITTj) with 

ds" = -^i-dt^ + {dx'f) + //(c/r^ + r^dnl) , (13) 

{i = 1, 2, 3) and 

^(5) = Q(e(5) + *e(5)), (14) 

where dfl^ is the metric on the unit 5-sphere, and e(5) is the volume 5-form on it. The 
metric function in (fT3|) is 

/ = 1 + ^ (15) 
and the relation between the various parameters is 

Q = lQ7rgs{a'yN = 4R\ (16) 



7 



where gs is the closed string couphng. The D3s themselves are hidden behind the 
degenerate horizon of f[T5]) at r = 0. 

Note that we are using the full D3 geometry, rather than the near- horizon limit; 
the latter would correspond to neglecting unity in comparison with R^/r"^ in (fT5|) . 
The physical significance of this has been discussed in the introduction. Technically, 
we need the full geometry because a probe D5 placed in it will be repelled to large 
values of the radius r, ~ R- 

The probe D5 is governed by the action 

Sd5 = -r5 / d'^^^J-detiPigU + Fa,) + r, J A A P[G(5)] , (17) 

where ts is the brane tension, with a = 0, . . . , 5 are coordinates on the brane, P 
denotes pullbacks to the brane worldvolume, Aa is the worldvolume gauge field, and 

Fab = daAt - dbAa (18) 

is its field strength. The WZ term is the second term in flTTl) . 

In our case, the D5 wraps only three of the five angles appearing in f[T^ . so it is 
convenient to rewrite the metric as 

ds^ = ^7 (-dt^ + {dx'Y) + Jf (dp^ + p^dnl + dA^ + A^dcj)^) , (19) 
vJ 

where A and are polar coordinates in the (x^,x^) plane. Then, 

= p2 + A^ (20) 

In our calculation, the instability with respect to the D5 slipping off the equator of 
the S"^, i.e., developing a non-zero expectation value of A, corresponds to a pairing 
instability towards a superconducting state. 

With these notations, the coordinates on the brane are 

r = (t,x,p,a), (21) 

where x = x^, and ex. = {a^, o? ^ a^) are the three angles spanning the in (fT9|) . If we 
restrict attention to embeddings with = = 0, the general form of the embedding 
is A = A(^"), = 0(^"), Aa = Aa{^^). It is, however, consistent to restrict the class 
of embeddings further, to 

A = A(p), (22) 
= gx, (23) 
At = A(p), (24) 
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where g is a constant, and all Aa with a ^ t are equal to zero. 
For this class of embeddings, the WZ term in ([T7|) becomes 

S^, = nJdtdxdpd^aAP[G,,,U^. (25) 

Only the magnetic part of f[T^ contributes to fl2Sl) . The requisite components of G'(5) 
are 



GA0c.(p,a,A) = -Ap'dpf{p,A)^gs{(x), (26) 
Gp^„(p,a,A) = Ap^dAf{p,A)J^), (27) 



G,A«(p,a,A) = 0, (28) 
where is the determinant of the metric on and 

/(P^^) = l+ (^2^1^)2 (29) 

[cf. (fTSl) and (l20l)]. The requisite component of the pullback is 

— g (Ga^qA p + Gp^a) , (30) 

where we use the shorthand A^p = dpA. Substituting this into ( l25l) and integrating 
over Q!, we obtain 

Swz = 27r2gr5i?' / dtdxdpAt{p)^Tl{p, A{p)) , (31) 
J dp 

where 

The full D5 action (H?]) for embeddings of the form (l22])-(l24!) is 

Sd5 = -2n\ J dtdxdpy/Z{l + A^ - F^^f^ + Swz , (33) 

where 

Z(p, A) = pV(p, A) [1 + g2AV(p, A)] , (34) 
with / given by ( p9|) . and A by A(p). 
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4 Numerical solutions 



We call an embedding solution superconducting if it is nontrivial (i.e., A(p) is not 
identically zero) and, in addition, 

A(p^oo) = 0. (35) 

The latter condition means that that there is no superconductivity in the ultraviolet! 
If the D5 brane does not cross the horizon, i.e., 

A(p = 0)^0, (36) 

the minimal distance from it to the D3s is finite and so then is the quasiparticle gap. 
We refer to such solutions as gapped. As we will see, there is a gapped supercon- 
ducting solution for any value of q. This is in contrast to the same system at fixed 
current [5], where the superconducting solution at any nonzero current is gaplessj^ 

The case g = is 2-dimensionally Lorentz invariant and so is subject to Coleman's 
theorem [15] on nonexistence of Goldstone bosons in 2 dimensions. This implies that 
the long-range order (LRO) represented by the nontrivial profile of A is destroyed by 
quantum fluctuations. One may suspect that the same is true also for the Lorentz non- 
invariant case q ^ 0^ The strength of quantum fluctuations, however, is suppressed 
by the large N (essentially, by the D5's tension in the gravity description or by the 
thickness of the wire). As a result, as increases, significant deviations from classical 
behavior occur at progressively larger spatial scales. Thus, our classical solutions 
faithfully represent the physics except at these largest scales. We also recall that 
superconductivity does not require an LRO; it only requires that the superconducting 
density does not renormalize to zero at large distances. 

The equation of motion obtained from (133|) by variation with respect to At inte- 
grates into 

"^^'^ + qR'U = - Jo , (37) 



{1 + Ai-Ftiy/^ 

where 11 is given by fl32|) . and Jq is the integration constant. Solving this algebraically 



for Ftp, we obtain 



Ft, = -^{l + Al)y\ (38) 



§In that respect, perhaps a better, if a bit cumbersome, name for these embeddings would be 
spontaneously superconducting. 

^The solution at zero current is gapped and coincides with the q = solution found here. 

II There is a subtlety, though. The gapless superconductor of [5] is described, at least near the 
critical current, by the dissipative XY model, which supports an LRO (at zero temperature) even 
in 2 spacetime dimensions. Qualitatively, that happens because low-energy fluctuations of the order 
parameter are damped in the presence of gapless quasiparticles. We assume that a finite gap prevents 
this mechanism of stabilizing an LRO from operating in the present case. 
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where 



J(p,A) = gi?^n(p,A) + Jo, (39) 
C(p,A) = Z(p,A) + j2(p,A). (40) 

The condition (!36|) requires that Jq = 0: unless that is so, Ftp{p = 0) 7^ 0, and At 
cannot be smooth at p = 0. 

Substituting (1351) into the equation obtained by variation of ( 1551) with respect to 
A, we obtain 

d_^/C^ = (1 + A^ )V2^ (41) 
dp(l + A2^)V2 + ' ^ ^ 

where C is given by ( HOl) (and J by ( l39l) . with Jq = 0). Regularity of A(p) at p = 
together with (136|) leads to the boundary condition 

A,p(p = 0) = 0. (42) 

Thus, (HT]) has to be solved with the boundary conditions (135|) and fH2|) . We do that 
numerically by shooting from p = with A(p = 0) as a shooting parameter. 

Some representative solutions are shown in Fig. [1] We find that such a solution 
exists for any g, no matter how large. The solution approaches a fixed shape as 
q — >■ 00. The absence of a classical stability bound on q (i.e., of the Landau criterion) 
has been remarked upon in the introduction and is further discussed in the conclusion. 
The minimal value of 

r{p) = \p' + /\\p)Y'\ (43) 

which gives the quasiparticle gap in units of the string tension, is reached at p = 0. 
Thus, the gap is given by A(p = 0). A somewhat counterintuitive result is that it 
grows with q (and approaches a constant value at g — )■ 00). 



5 Worldvolume flux and charge 

The counterpart of (!37|) for the fixed current case [5j is obtained by setting g = 
and retaining a nonzero Jq (which is then interpreted as the fixed current). Let us 
compare the fluxes of Ftp at p — )■ 00 for the two cases. They are given by 

Ftp{p^ 00) = (fixed g), (44) 

Ftp{p 00) = (fixed current). (45) 
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Figure 1: D5 profiles for tliree values of q. Both axes are in units of i?, and q is in 
units of 1/R. Larger q correspond to larger values of A(p = 0). 

These are obviously similar. Indeed, as we will see shortly, qR^ can be interpreted 
as the average current in the fixed-g case. Note, however, that, unlike ( H5l) . which is 
sourced by charges behind the horizon, the flux ( H4l) is entirely due to charges on the 
brane. 

According to fl3T|) . the total charge, coupling to the properly normalized potential 

At/(2W), is 

An'^a'T^R'^ J qdx = NW , (46) 

where W = qL/2TT is the winding number, and we have used the expression f lTB]) for 
R and [I6] 

= {2nfgs{a'r ^^^^ 

for the brane tension. Note that (146|) is an integer. In the context of the original 
condensed-matter system, we identify it with the total momentum carried by the 
current, in units of the Fermi momentum. 

Consider a process (a quantum phase slip) in which the winding number W de- 
creases by unity. According to fH^ . this releases N units of charge from the brane. 
Since the charge is conserved, some other charged objects must appear. For a differ- 
ent brane arrangement, where a D3 and a D5 have no common spatial dimensions, it 
is known that a fundamental string stretching between the branes is produced when 
they cross; this is the Hanany- Witten effect [S] in one of its dual versions [HI [IDl E] ■ 
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Charge conservation leads us to conclude that the same is true here. In the present 
context, a fundamental string corresponds to a quasiparticle in the superconductor, 
so we conclude that each phase slip must produce N such quasiparticles. This is 
precisely the result we have obtained previously by instanton computations on the 
basis of the Bogoliubov-de Gennes equations [7]. 



6 Computation of free energies 

One consequence of quasiparticle production by a phase slip is that it affects the 
energy balance: the energy gain from unwinding the supercurrent must be weighed 
against the energy cost of the produced quasiparticles. As a result, when production 
of quasiparticles is a requirement, phase slips may be energetically forbidden (at zero 
temperature), and the supercurrent stable. In this section, we describe a numerical 
computation of the relevant free energies. 

Substituting (155]) into the expression ( 155]) for the action and integrating by parts 
in the WZ term, we obtain 



= 27rV5 J dtdx - J dpVC{l + Ay + qR'^Atioo) 



(4J 



The integral over p here is divergent at p -^■ oo, but that can be fixed by subtracting 
some reference independent expression, e.g., 

VC{1 + Ay 1/2 ^ Veil + A2^)1/2 - y^;, (49) 

where Cq = p^(l + i?Vp^). 

The asymptotic value At{oo) appearing in (1^8]) can be interpreted as a chemical 
potential for the electric current. This can be seen as follows. The average current 
carried by electrons near the Fermi surface is given by their momentum density in 
units oi kp. In the D-brane description, the momentum density is represented by 
the density of charge conjugate to At. Thus, changing the chemical potential for the 
current corresponds to a shift 

A(p)-> A(p)+Ao, (50) 

where Aq is a constant. Allowing the D5 to have a nonzero total change implies that 
the physical states of the D5 are only invariant under gauge transformations with 
parameters vanishing at infinity, in particular, At{p) — )■ At{p) + A(p) with A(p) — )■ 
at p — 7- oo. Then, At{oo) is a gauge-invariant quantity (a collective coordinate in the 
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terminology of [E]), precisely the one mo dified by the shift ([50]) Q 

Variation of (HSj) with respect to At{oo) gives the average current / = 2t{'^t^J, 
where 

J = qR\ (51) 

Legendre transforming the action with respect to At{oo) and peeling off (— / dt) gives 
the free energy of the supercurrent state 

^super(g) = 27T%Le{q) , (52) 

where 

Q{q) = / dp 
Jo 

We now consider various final states to which supercurrent may possibly decay. 
The first group includes two gapless states: the normal state A = and the gapless 
superconductor found in [5j . A decay to either of these can be visualized as production 
of a finite density of strings that pull the D5 through the horizon, hiding behind it 
the worldvolume charge that has originally resided on the brane. The momentum 
(charge) conservation implies that, to see if either decay is energetically allowed, we 
must compare the free energy (152!) to that of the final state with the same value of 
the electric fiux at infinity. According to ( 14^ and ( H5|) . this corresponds to setting 
Jo for either state equal to J, eq. ( 151]) . of the supercurrent state. The free energies of 
the gapless states are subtracted in the same way as (155]) : thus, for instance, the free 
energy of the normal state is 

J^novmiq) = 27rV5L / dp 

Jo 

Numerically, we find that either of these free energies is always higher than that of 
the supercurrent state, so the decay is energetically forbidden. 

Another type of final state is a gapped supercurrent state with a smaller value 
of q. Here we limit ourselves to states connected by a single phase slip, AW = — 1, 
in a wire of macroscopic length L. In this phase slip releases approximately 

{2tt / L)dJ^supcr/ dq of free energy from the supercurrent. It, however, requires pro- 
duction of N quasiparticles. The minimal energy of a quasiparticle is given by the 
quasiparticle gap, (27ra')~^A(p = 0). It is convenient to adopt the system of units 
where all lengths are measured in units of R, i.e., R = 1. In these units, the co- 
efficient in (15 2 p can be written as 27r^r5 = N/An'^a'. So, to see if a phase slip is 

**Note that, unlike other identifications of chemical potentials in D-brane systems (for example, 
the one in [18] ) , the present argument does not rely on the AdS / CFT dictionary of [11 |3l |4] . That 
dictionary cannot be used here as p — oo corresponds to the asymptotically flat region of the 
geometry. 




(Co + Ji 
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Figure 2: The amount of free energy (in units of NR/27ra'; solid line) released by a 
single phase slip from the supercurrent in a long wire, as compared to the quasiparticle 
gap (in units of R/27ia'; dashed line), both as functions of the winding number density 
q. Intersection of the curves marks the value of q (in units of 1/-R) at which a phase 
slip becomes energetically allowed. 



energetically allowed, we need to compare dQ/dq to A(p = 0). The comparison is 
shown in Fig. O We see that the decay becomes possible for q larger than a certain 
qm (which is somewhat below 0.6, in units of 1/-R). 



7 Discussion 

Overall, we find that supergravity provides a remarkably detailed picture of a clean 
(disorder-free) multichannel one- dimensional superconductor. The picture includes 
the requirement of quasiparticle production by quantum phase slips, seen here as a 
version of the Hanany-Witten effect (in a non-transverse, non-supersymmetric ar- 
rangement of branes). This complements the earlier derivation [7j based on the 
Bogohubov-de Gennes equations. 

Gapped classical solutions exist in the present case for any winding number density 
q, no matter how large. As we have noted in the introduction, this may be a conse- 
quence of momentum conservation, which precludes supplying a nonzero momentum 
to the superconductor as a whole (the Landau process). It would be interesting to 
see if there is a limit to classical stability in a case when momentum conservation is 
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broken. One such case occurs when the stack of D3s is interrupted by an NS5 brane, 
in an arrangement similar to those considered in The numbers of D3s on the two 
sides of the NS5 need not be equal, so placing a probe D5 in such a geometry will 
represent a wire with a varying number of transverse channels, i.e., a constriction. 
The Landau process in this scenario would correspond to formation, at a critical cur- 
rent, of a gapless region near the NS5. It remains to see, of course, if that is indeed 
what happens. 

Unlike the Landau process, a quantum phase slip does not change the total mo- 
mentum of the system: it only transfers momentum between the supercurrent and 
quasiparticles. We have seen that the supercurrent-only solution is metastable beyond 
a certain g = g^. If we try to set up in a ring a supercurrent larger than = QmR'^, 
a part of it will decay into quasiparticles and, once these form a Fermi surface, a gap- 
less normal component will appear. The latter will exist alongside a superconducting 
component — as we have seen in Sec. [6l in the present perfectly uniform wire 

at zero temperature), the normal-only state is never the most energetically favorable. 

In the presence of disorder, quantum phase slips can occur without quasiparticle 
production Nevertheless, our present results lead us to consider the possibility 
that, in that case too, a sufficiently large current (now maintained by an external 
battery) causes appearance of a gapless normal component in the superconducting 
state. One may contemplate trying to detect such a component experimentally — for 
instance, by measuring the current-voltage curve of electrons tunneling into the wire 
off the tip of a scanning tunneling microscope. 

This work was supported in part by the U.S. Department of Energy grant DE- 
SC0007884. 
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